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We explore the effect of two-dimensional position-space non-commutativity on the bipartite en- 
tanglement of continuous variable systems. We first extend the standard symplectic framework for 
studying entanglement of Gaussian states of commutative systems to the case of non-commutative 
systems residing in two dimensions. Using the positive partial transpose criterion for separability 
of bipartite states we derive a new condition on the separability of a non-commutative system that 
is dependent on the noncommutative parameter 9. We then consider the specific example of a 
bipartite Gaussian state and show the quantitative reduction of entanglement originating from non- 
commutative dynamics. We show that such a reduction of entanglement for a non-commutative 
system arising from the modification of the variances of the phase space variables (uncertainty 
relations) is clearly manifested between two particles that are separated by small distances. 
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I. INTRODUCTION 



It is a common expectation that spacetime may not retain its smooth continuous structure at very short distances, 
and hence can have important consequences for physics at very high energies, like quantum gravity. Moyal deformation 
of ordinary spacetime is an example of a specific algebraic realization of this expectation [l[ . Clarification of the role of 
underlying spacetime symmetries is crucial for formulating quantum dynamics on such spaces 0, &@|- However, 
it is not only for considerations of quantum gravity that the study of non-commutative spacetime is relevant: non- 
commutative coordinates also arise in low-energy situations as well, like in the Quantum Hall Effect 01 • In the 
presence of a magnetic field, the guiding centre coordinates of the electron become non-commutative. Understanding 
the formulation of quantum theory on such a space is irnportant to address basic conceptual issues like the behavior 
of many-particle systems, notion of identity of particles 0, and connection between spin and statistics. 

While it is not entirely surprising that spacetime non-commutativity may lead to non-locality at short distances, 
we wish to address another novel aspect of non-commutative quantum mechanics here. How are notions of entan- 
glement and information affected ? Quantum entanglement is the key ingredient for the storage and distribution of 
quantum information among the fundamental constituents of the quantum world @ . The first nontrivial consequence 
of entanglement on quantum ontology was noticed many years ago within the context of position-momentum con- 
tinuous variable systems Q. In recent times there has been a rapid development of the theory of entanglement for 
such systems Several interesting information processing applications such as quantum teleportationlll| . dense- 
coding [l^ and cryptography 13 1 have been extended to the domain of continuous variable entangled states[l4 fisj. 

Understanding the physical origin of entanglement in quantum systems has been a key goal of information theory, 
which is expected to lead to insights on a vast range of diverse phenomena such as phase transitions in condensed 
matter systems [l^ and black hole physics [13 • Such a goal is unlikely to be completely realized unless the role of 
the symplectic structure of the underlying phase space on the information capacity or entropy of quantum systems 
is clearly elicited. It seems likely that any modification to this structure, as brought about by position-space non- 
commutativity at a fundamental scale could impact the correlations of quantum systems. There indeed should exist a 
deep-rooted connection between non-commutative dynamics and the phenomenon of entanglement of quantum states. 
Our motivation for this work is to initiate a systematic investigation on the role of position-space noncommutativity 
on entanglement. In the present st udy we restrict ourself to the entanglement of bipartite Gaussian states for which 



tudy we res 

the standard symplectic formalism [Tsl. [Tol [20| is readily available and well-established in the commutative limit. 
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The plan of this paper is as follows. In Section IIA we briefly review the standard symplectic formalism for 
determining the separability of bipartite Gaussian states in a single spatial dimension. In Section IIB we first motivate 
the need for considering dynamics on a two-dimensional plane relevant for non-commutative systems and then show 
how the variances are modified under the effect of non-commutativity leading to ^-dependent uncertainty relations. 
In Section IIC we extend the usual commutative one-dimensional (spatial) formalism for application of the positive 
partial transpose [Ull^l separability criterion for Gaussian states [13 to the case of a bipartite system in two (spatial) 
dimensions. Using the symplectic formalism we derive the modified (0-dcpendcnt) condition for the physicality of 
states on a non-commutative plane which leads to a corresponding modification of the entanglement condition. We 
then consider various examples of Gaussian states. In Section IIIA the entanglement of a commutative bipartite system 
in one (spatial) dimension is computed. In Section IIIB a single particle on the two-dimensional noncommutative 
plane is studied showing that the underlying non-commutativity does not impact the a priori separability of the 
modes in this case. In Section IIIC an example of a bipartite state in two dimensions obeying the usual commutative 
dynamics is considered in order to illustrate the validity of the symplectic formalism in two spatial dimensions. We 
finally investigate a bipartite Gaussian state on a two-dimensional non-commutative plane in Section HID. We show 
that the 0-dependendent entanglement criterion leads to a reduction of the magnitude of entanglement. We conclude 
with a summary of our results in Section IV. 



II. EFFECT OF NON-COMMUTATIVE DYNAMICS ON ENTANGLEMENT OF GAUSSIAN STATES: 

FORMALISM 



A. Separability condition for a bipartite system 

To be gin with, we will first review the symplectic formalism for observing entanglement of a bipartite Gaussian 
state[l8lll9l[20| where the usual position space commutativity is maintained. This analysis is in order not only to make 
the paper self-contained, but will be useful as a ready reference when we extend it to incorporate non-commutativity 
in a subsequent section. Let us consider a composite wave function ^' describing the bipartite state of two particles 
(in momentum space) as 

^{pi,P2) = Ni[^Pi{pi) ® ^P2iP2) + MPi) MP2)] (1) 

where A'^i is a normalization constant. 

We have to first construct the variance matrix for the composite wave function 5*. We can then apply the positive 
partial transpose (PPT) separability criterion. The PPT criterion is a necessary condition for separability of quantum 



states of arbitrary dimensions [2l|. It is also a sufficient condition for bipartite 2® 2 and 2 (g) 3 systems. Here we use 
the generalization of the PPT criterion for Gaussian states [l^, [l^] based on the smallest symplectic eigenvalue of the 
partially transposed variance matrix. The variance matrix elements are calculated from the variances and covariances 
of the phase-space operators xi,pi, X2,P2 hi the state ^'(p). The ordering of the phase-space variables in the variance 
matrix for the state ^'(p) is chosen to be = {xi,pi,X2,P2)- The 4x4 matrix formed by the set of bilinears C^<^^ 
transforms under symplectic transformation as 

} - } = {Si^S^ti''} = S{ti^}S^, S e Sp{A, R) (2) 

where the indices /i, i/, a, (3 range from 1 to 4. A general element of the operator matrix can be written in terms of 
the commutator and anti-commutator as 

= + (3) 

where the real variance matrix V for the state p is defined by 

m = Tr(pee^) 

^ V + i^l (4) 

where 

y = K''^iTr(p{e^r})} (5) 
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can be identified as the variance matrix and the matrix 



(v Trp = 1) (6) 
is the symplectic matrix. For a single mode system the variance matrix is given by 

For any matrix to be physicaUy reaUzable (respecting the uncertainty condition) one requires that 

dotV ^ {x^)lf) - [\{{xMf >\ (8) 

(in units of h ~ 1). This is the symplectic invariant form of the Heisenberg uncertainty relation AxAp > 1/2, 
where we have assumed without loss of generality that (i) = = (p) by suitable translation in phase space, so 
that (z^) = (Ax)^ and (p^) = (Ap)^. Otherwise, the entries of the variance matrix will take the form V = {v^'^ = 
^((C^C ^ which are clearly invariant under phase space translations. The variance matrix V can 

be diagonalizcd by suitable symplectic transformations to the form V ~ diag{i' /2,v /2\ . The physicality condition 
then takes the form 

> 1 (9) 

For an n-mode system, let us consider the state /3(s) = U (s)pU^s), where J7 is a unitary operator implementing 
the symplectic transformation S £ Sp{2n,R). It then follows that [1^ 

V' + svs'^ + in (lo) 

It is clear that if y is a physically realizable state, so is its symplectic transformation V' for every S € Sp{2n,R). 
Further, using Williamson's theorem [23|. it is possible to transform the variance matrix V using symplectic transfor- 
mations to the canonical scaled diagonal form, up to the ordering of i^j , given by 

V' = SVS^ - diag{iyj2, 1/2/2, • • • , z/„/2, i/i/2, 1^2/2, • • • , i/„/2) (11) 

This set of eigenvalues {lyj} which are at least doubly degenerate, forms the symplectic spectrum of V and can be 
obtained alternatively by computing the ordinary eigenvalues of |2ir2V^|. One can then treat the n-modc system simply 
as n copies of the one-mode system, whereby V and V are bona fide (respecting the uncertainty relation) variance 
matrices, if and only if (iff) 

i^]>l, .7 = l,2,....n (12) 

The matrix given by V + in has a spectrum of eigenvalues (usual eigenvalues, not symplectic eigenvalues) {vj ± l)/2 
for j = 1, 2...n. It follows that the physicality condition (|12p amounts to demanding that V + in> 0. The symplectic 
transformation 

V + in = s{v + in)s^ (is) 

connecting the two matrices {V 4- ifl) and (V' + ifl) preserves the condition of positive scmi-definitencss, i.e., 

V + in>0 (14) 

In order to use the PPT criteria (as applied for continuous variable systemsfl^ [2^) for a two-particle state, 
the variance matrix V is transformed into V by applying time-reversal symmetry on the second particle (partial 
transposition). Equivalently, this is obtained by flipping the sign of the momentum of the second particle so that the 
variance matrix V undergoes the following transformation 

V = AVA (15) 

where A = diag{l, 1, 1, —1). We then again find the eigenvalues i^i of the matrix |2ir2F|, which form the symplectic 

spectrum of the variance matrix V, where {I = (B uj, oj = For separability of the state, the PPT criterion 

reduces to a simple inequality satisfied by the smallest symplectic eigenvalue i>_ of the partially transposed state [l^ 

i>l > 1 (16) 

Any two- mode Gaussian state is separable if and only if the above inequality holds, otherwise it is entangled. 
The magnitude of entanglement can be estimated through the logarithmic negativitvpsj defined as 

E = inax[0,~\og^{D_)] (17) 
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B. Kinematics on a non-commutative plane 

The algebra satisfied by the operator- valued coordinates in the simplest form of non-commutative space is given 

by 

=i6i^^ (18) 

where 9^^ is usually taken to be the non-transforming entries of the constant antisymmetric matrix Q ~ {0'"^}. 
In the usual 3 + 1 dimensional spacetime, the introduction of invariant length scales through Q, however creates 
problem regarding Lorentz symmetry, which can be restored only in a twisted framework in a Hopf-algebraic setting 
[2]. However, as it can be easily seen that for the special case of a 2 -|- 1 dimensional spacetime, with time t taken to 
be a commuting variable. 

[xi,X2]=i0, = [£,l2] =0 (19) 

The entire set can be written covariantly (under S0{2)) as 

[xi,Xj\ ^ ideij ^ iOij, [£,lj]==0 for 1,2 (20) 

where 6ij = 6eij, {cij is the Levi-Civita alternating tensor with ei2 = -fl) representing the spatial parts of the matrix 
0, transforms as a scalar under SO{2). The space part of this 2-f 1 dimensional s pac etime will henceforth be referred 
as non-commutative plane, for which 9 is the single non-commutative parameter [23|. 
Our basic commutation relation therefore is 

= iOij = iOcij. (21) 

Instead of working with operators ij, one can work with c- number valued coordinates Xi, provided one replaces 
operator products with Moyal star multiplication, which is defined for a pair of generic functions f{xi) and g{xi) 
belonging to Schwarz class [2y| 

1 _ - 9 

f{x)-kg{x) = f{x)eyip{-e€ijdid,)g(x), d, = — (22) 

2 oxi 

This can be expressed alternatively as 

f{x) * g{x) = mg{f{x) (g> g{x)) = mo(Fe/(S) «) gix)) (23) 

where Fg = exp(^9eijPi®pj) is the Drinfeld twist operator and mo represents pointwise multiplication i.e. mo(/(x)(8) 
^(a;)) = f{x)g{x). With this, the Moyal bracket 

becomes isomorphic to Eq. (PT|) . This generates the non-commutative algebra y^e(M^) and represents a deformation of 
the corresponding commutative algebra ylo(]^^) where the functions are composed through point- wise multiplication. 
Because of this commutation relation, making boxes in noncommutative space is tricky at best, so we will avoid this 
route. Instead, we work directly with wave packets, characterized by its widths and average momenta/locations. But 
before doing so, we must first understand how to do quantum mechanics on the non-commutative plane. 

In usual quantum theory, the i*'' coordinate position is represented by an operator Xii it acts on a wavcfunction 
ipix) in the position representation by multiplying ip{x) by xf. 

Xi'ip(x) = Xiip{x). (25) 

For the non-commutative plane, we need to distinguish between left - and right-multiplication operators: 

xf-ip^x) = Xi-kil)(x), (26) 
-0(5;) = il}{x)-kXi. (27) 

For concreteness, we shall henceforth choose x^ as operators corresponding to coordinate observables Xi. 
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Our basic quantum mechanical commutation relations are: 



It is easy to check that 



i^i, Pj] = i^i]! 

\p,,p,] = 0=[^f,i^]. (28) 



1 ^R. , . 

Pk = -^f-k,n{x,„, - X„^) (29) 



satisfies the above commutation relations. In fact, puipix) — —idkip{x). Note that this is the same as the adjoint 
action of momentum introduced at the level of operators in . 
Notice that 

^ i±i (30) 

form a commutative quantum algebra: [xi, xj] = 0. So the quantum operators Xi can be simultaneously diagonalized, 
and their simultaneous eigenvalues generate the commutative algebra ^o(]^^)- 

Also, one can easily check that [xi,pj] = iSij, so Xi and pj generate the usual Heisenberg-Weyl (HW) algebra. 

Xiip{x) — Xiip{x) (no star!), (31) 

1 -L -R — 

Pi'^ix) = -e^jn{x^-x^)'4!{x) = -idiip{x). (32) 

I ■ ■ ^ = "} = {xi,X2,pi,P2)'^ be the corresponding phase-space coordinates. The complete NC Heisenberg 
algebra is obtained by augmenting Eq. (j2ip by the following additional commutators 

[xi,pj] = i6ij, [Pi,P]] = (33) 

Let the corresponding underlying commutative structure be described by the phase-space coordinates ~^ = {^'^} = 
{xi,X2,pi,P2)^ which satisfies the usual (commutative) Heisenberg algebra, i.e., for which 6* = in Eq. ipTj) . These 
phase-space coordinates and ^" are related by the following linear transformation 

e = M'' X (34) 

where the matrix M is given by 



(35) 



Clearly, their NC counter-parts V and Q. are related to V and Q. by 

V = MVM^ and Vl = MTIM^ (36) 
Let us now consider the following wave packet in momentum space: 




1 



V'(p; a,Po, a) = -^eS-^^P-^^ -'P'^^ (37) 



where bi = ai -\ — ^iM^^h. ^ ^j^g reason for this particular choice of the form of b will become clear subsequently. The 
representation of Xi in the momentum space is given by 

Xii'{p;a,po,a) = {idp, - ^^^)V'(p; a,Po, «) (38) 
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We would now like to calculate the spreads in Axi, Ax2, Api, Ap2 in the state given by Eq. 
arguments one expects that 



37|) . From general 



(Aii)(Ax2) > 



(39) 



(ASi)(Api) > - 



(40) 



(Al2)(Ap2) > 2 



(41) 



To derive the uncertainty relations, we proceed with the variance matrix approach. The elements of the variance 
matrix in the basis of the phase-space variables ^ = {ii, X2,Pi,P2} in the non-commutative plane is obtained as (27j 



where 



^ Mr 
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4 
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2 / 



(42) 



(43) 



For example, it turns out that (x)^, = a, (p)^ = indicating that the particle represents by packet (j37p is located 

at a and has a mean momentum ^ . These simple forms for the mean position and momentum explains the choice of 

the particular form of h in Eq. ([37|) . 

From the above variance matrix, we have 



Axi = Aa;2 



612 a 



1 



Api = Ap2 



8 2a 

Therefore, the "space-space" uncertainty in the state ip{p;a,po,a) is given by 



This uncertainty is minimized when a = 



Aa:iAx2 = —— -I- — 
8 2a 

: I, and 



(44) 



(45) 



(AiiAi2)r^ 



(46) 



It follows that the minimum space-space uncertainty cannot be arbitrarily reduced below the value given by Ea. (|46p . 
This is different from the case of the standard commutative dynamics {9 = 0) where by increasing a (the spread of 
the wave packet in momentum space), one can measure the observables Xi,X2 to arbitrary precision. Similarly, the 
"phase-space" uncertainties in the state ipip; d,po, a) is given by 



AxiApi = AX2A.P2 = 



16 



(47) 



For a = ao, the "phase-space" uncertainties reduce to AiiAp2 = Ax2Ap2 = The above equations (|i5l) and 
P7)) determine the space-space and phase-space uncertainties under the effect of non-commutative dynamics. For 
instance, the variances are clearly enhanced in the non-commutative case {9 7^ 0). This suggests that by increasing 
the uncertainty, i.e., the fuzziness of the unit cells in phase space, the effect of non-commutativity is to decrease the 
available degrees of freedom (the total number of unit cells) if the overall volume of the phase space is held fixed. 
The number of allowed states in a noncommutative system is essentially ^ ^ for a planar system of large but finite 
area A. The fact that this number is reduced in presence of non-commutativity, was utilized to address the issue of 
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fractional Quantum Hall Effect in Ref. [23| in a fermionic s^tem involving electrons, where the entropy gets reduced 
generically, which also displays non-extensive behaviour [29[. One therefore expects a similar reduction in entropy 
in the bosonic case also, just on the basis of enhanced phase-space uncertainity (|47p. which in turn is expected to 
play a vital role in the entanglement property as well. The consequences on the entanglement of two particles will 
be evident later when we present explicit computations of the entanglement of bipartite systems in the presence of 
non-commutative dynamics. 

Note that the commutative 6* — > limit of the variance matrix V^^j reduces to Vq = diag{^, ^7 f 7 f )■ The same 
form of this commutative variance matrix Vq can also be obtained alternatively from by making use of the linear 
transformation (j36p to express it in terms of the commutative coordinates. The various uncertainty relations can now 
be read off from this as 

AxiAx2 = ^, AiiApi = A.T2AP2 = J (48) 
Za 2 

Clearly, by taking the spread a of the momentum space Gaussian wave packets arbitrarily large a ^ 00, we can find 
(Ail Ax2)min = 0. This disappearance of the 9 — dependence from Vq, or the various uncertainty relations following 
from it, is a reflection of the fact that one does not expect to observe any non-commutative effect at the single 
particle level, if the commuting variables are used. We shall see subsequently that this feature does not persist in the 
presence of two or more particles, as the symplectic eigenvalues of the effective commutative variance matrix will start 
displaying 0-dependence (non-commutative effect), thus inheriting the original non-commutative features. However, 
in the limit of a — s- cx3, the space-space uncertainty relations involving non-commutative coordinates diverges (j45p : 
(Aa;iAa;2) oo, which is minimized only at ao = | as we have seen earlier. 



C. Entanglement of two particles on a non-commutative plane 



For a commutative system the correlation properties of the two-dimensional Gaussian state is completely determined 
by an 8 X 8 real symmetric variance matrix V . To study the correlation properties it is convenient to transform the 
variance matrix V to some standard form through symplectic transformations. Let us choose the basis vector for the 
phase space variables to be ^ = (x^^^p^^'', ij^^pj^^ x^^^pp'', a;2^\p2^'), where the subscripts 1 and 2 refer to the x 
and y components, respectively, and the superscripts 1 and 2 (within parenthesis) refer to the first and second particle, 
respectively. Any variance matrix can be written in terms of 4 x 4 block matrices a, (3 and 7 as 



a 



8x8 



where 
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(49) 



(50) 



represent the single particle (co-) variances for the first and second particle, respectively, and 
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represents the correlations between the two particles 
We now extend the procedure presented in Rcf. [20 
any general variance matrix ()49p to a standard symplectic invariant form 
procedure). The standard form of the variance matrix is given by 



(51) 



to the case of a (spatial) two-dimensional system to transform 
(Refer to Appendix I for details of this 
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(52) 



8x8 



8 



A Gaussian state described by the variance matrix V is physical if the smallest symplectic eigenvalue ^_ > 1 with 
the constraint F > [H, Eol. 

The stage is now set to formulate our prescription for studying the entanglement of a non-commutative system. 
To discuss the entanglement properties of two identical particles in two-dimensional non-commutative space, we must 
first define such states appropriately. We briefly recall here the arguments provided in [1, Q that give the correct 
construction. 

In noncommutativc spacctimc, generators of symmetry transformations act differently, because this action has to 
be compatible with the star-product. The correct implementation of, say, the rotation operator J requires that J act 
on a two-particle state not via the usual product Ao( J) = 1(E) J + J but through the deformed co-product defined 
as, Ag{J) = F^^Ao{J)F0, where Fg'^ = exp[(^)(pj (X)_Pj)] is the inverse of the twist operator. 

In usual quantum mechanics of identical particles, all observables commute with the statistics (or flip) operator 
To that exchanges the two particles: tq is supers elected. The twisted two-particle rotation operator A.g{J) does not 
commute with fo, but does so with the twisted flip operator fg = Fg^^foFg = F^^'^to (An easy way to verify the second 
equality is to look at the action of fg on, say, plane waves). 

To ensure that statistics remains superselected, we therefore require that the two-particle wave function be 

* = 7V(?Ai®V2 + -Fe"'(V'2®^i)) (53) 

where N is the normalization constant. Such a wavcfunction transforms correctly under the action of spacetime 
symmetries (roughly speaking, the wavcfunction remains twist symmetric in all reference frames). 

It has been shown recently that these twisted particles can have some non-trivial consequences like violation of 
Pauli principle [sO] and energy shift in the system of degenerate electron gas [3l|. Though 6*^ = 9eij turns out to be 
an SO(2) scalar as mentioned earlier and the technique of twist symmetrization is motivated by its generalizability to 
higher dimensions. Indeed, general arguments coming from diffeomorphism invariance [s^ and quantum integrability 
[33| support this. 

How does one generalize the phase space operators to the two-particle sector? In standard quantum mechanics, the 
position operator in the two-particle sector is Ao(x) = xE)i + i(Ex, of which we intuitively identify the first piece 
X (g) 1 as the " position of the first particle" . This identification is strengthened by the observation that on applying 
the flip operator tq to this, we get 1 ® x, which is the "position of the second particle" . In the noncommutativc case, 

the phase space operators are the exact analogues of these: they are given by fg~^(i 8) C )Pe and Fg^{£^ ^ i)Fg, as 
one can be obtained from the other by applying fg. 

For notational convenience, one can formally divide the (8 x 8) matrix into four blocks, i.e., 

y _ { (^1 ® ^)4x4 {mixed terms)4x4 
y{mixed terms)4x4 {lE)Vi)4x4 

The left-hand upper corner block and the right hand lower corner block contains the elements oi Vi® I and I (EVi, 
respectively. The other two symmetric blocks contains the mixed terms. The above matrix V is obtained from the 
symmetric part of the matrix {F^^£^£_'^ Fg) , where = {xi® I,X2® I,pi® I,P2® 1,1 ®xi,I ®X2jI ^Pi^I ®P2)- 

The important condition of positive semi-definiteness that is preserved under symplectic transformations for com- 
mutative systems ()13|) . however does not generally hold true for non-commutative systems. This is because VL is 
preserved only by the matrix S = Mj"^ 571^2, i.e., 

Tl^STls'^ (55) 

where M2 is given by 

(M®I \ 

{ i®m) 

with is defined in Eg. psp (33 |. Thus, the strategy employed here will be to transform the variance matrix V, 
to the effective commutative variance matrix V using Eq. (j36p . which also corresponds to the (co-) variances in terms 
of observables, as discussed earlier in Section IIB. The rest of the analysis parallels the one used for commutative 
systems, discussed in Section IIA. As a result of the transformation the non-commutative parameter enters into 
the effective commutative variance matrix elements. Our prescription will be to first transform a general form of the 
noncommutativc variance matrix to an effective commutative variance matrix, and then use the standard form for 
the latter to compute its symplectic eigenvalues (see Eas. ([55|) - ([W)) below). In general, any non-commutative variance 
matrix of two particles when transformed to the effective commutative plane will acquire 0-dependence. 
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It follows from Williamson's theorem [2J] that the effective commutative variance matrix V can be transformed by 
symplectic transformations to the canonical scaled diagonal form where one can separate out the ^-dependence of the 
various terms. Thus, one can write 

v = svs^ = siy{B = ^) + w(e))s^ 

= diag{vi{e = 0) + 5vi{e), • • • , 1/4(6* = 0) + 5vi{e), vi{e = 0) + 5vi{e), • • • , 1/4(6! = 0) + 5vi{e)) (57) 

In the commutative limit {6 = 0), the uncertainty relation imposes the condition v^{0 ~ Q) > 1 on the eigenvalues 
(see, section IIA). It thus follows, that the effective variance matrix V is physical if and only if {vj{0 — Q) + 5vj{9)Y > 
1 + {Svj{0))^ + 2vj{e = Q).5v.j{e) Now, consider the matrix V + iVl. It can be shown that this matrix has a spectrum 
of eigenvalues given by [yj{9 = 0) + 5vj[9) ± l]/2, for j = 1, 2, 3, 4. Since the matrices V + iVL and V + are related 
through a symplectic transformation, it follows that the positivity condition is preserved, i.e., 

V + in>{) (58) 

if 5i'j{9) > which turns out to be true in the example we are considering here. 

We therefore transform the non-commutative variance matrix into the effective commutative variance matrix using 
the relations (jSG]) . We then proceed in the same way as for the commutative system. Due to the 6'-dependence of 
the symplectic eigenvalues of the effective commutative variance matrix, the condition analogous to the physicality 
condition (|12p for the commutative case is modified to a ^-dependent relation here, given by 

{u^.^f > (.5^(0))L„ (59) 

where {v^'~^ {9))rnin is the lowest symplectic eigenvalue of the effective commutative variance matrix when the space- 
space uncertainty /S.xi/S.X2{9) is minimized. Note that {v^'~^)m.in ~ 1 when 9 ~Q (since AiiAi2 = in this case), i.e., 
for the commutative case, one recovers back the usual physicality condition for commutative systems. As a result the 
similar condition on the smallest symplectic eigenvalue of the partially transposed variance matrix (PPT criterion) 
also gets modified, showing that the entanglement criterion for noncommutative systems 

also acquires a ^-dependence We will derive the explicit form for [v'^'^ {9))min in the context of the example that we 
will study in Section HID. 



III. COMPUTATION OF ENTANGLEMENT: SOME EXAMPLES OF GAUSSIAN STATES 



In this section we will provide examples of specific Gaussian states for which bipartite entanglement will be computed 
using the formalism presented in the previous section. Though our ultimate goal is to study the effect of non- 
commutativity on the entanglement of two particles, we will proceed towards this direction step-by-step. Hence, 
following the sequence of the previous section, we will proceed by considering examples of entanglement in usual 
commutative dynamics in one and two dimensions, before launching into the non-commutative computations. This 
is in order to illustrate systematically the formalism provided in the previous section. 



A. Two particles in one dimension obeying usual commutative dynamics 



Let us consider two particles in one dimension described by the two normalized wave functions Tpi and V'2 given 
respectively, by 

/ ( \ 1 , ~(p-po/2)% 

V'i(p) = -Tcxp( ), (61) 

I ( \ 1 , -(p + Po/2)% 

W) = — — rcxp — ) (62) 

(tt/?)! 2P 

where a and (3 denote respectively the spread of the two wave packets. If we now construct a joint two-particle 

state given by Eq.((T]), with Nl = [2(1 + -jj-i j-j-^g normalization constant, one can observe that this 

^(^11^2) two-particle state cannot be written as a product of two wave functions of the single particle states unless 
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the parameters a and /3 are equal and po = 0. Therefore the state ^'(j)i,p2) is entangled when a 7^ /? or when a = P 
and Pq 0. 

Explicitly the variance matrix for two particles in the state \['(pi,p2) in Eq.([T]) is obtained as 



V 



/ Vii V12 Ui3 

W12 V22 V23 ^^24 1 /■ a' 7' 

wi3 W23 W33 U34 ) ~ V7'^ P' 

Vl4 V24 V34 V44 



(63) 



where a', /?', 7' are 2x2 block matrices. Using the tensor product notation for = (f (g) /, p (g) /, / (g) a;, / (8) p), one 
gets 

'a + B 4^/;^^exp(^) p2 



W12 = W34 = ((a; P + p i)/2 ® /)* = (/ (g) (a; p + p i;)/2)* = 



Wi3 = (g i)* = A^f 



4p2a/?cxp(^) 



/^(a + /3)3 

wi4 = W23 = (i g)p)^ = (p g) i)* = 

/^2^n ,.^^2. ,,2 + , , 4^/^exp(g|) gp pl{f3 - a)^ 

V22 = V44 = (p g) /)* = (/ g)p = — ^ ^ 3 \ 



a + P a + P 4:{a + py 



^2 pl^(p-afeM^) 



^;24 = (p®p)*=iViM + ^^^p ^1 (64) 



Needless to say, the above variance matrix (|63p corresponds to a physical state, and it can be easily verified that the 
eigenvalues of V satisfy the physicality condition Uj > 1. 

In order to study the entanglement of the state described by the variance matrix V, we now perform a partial 
transposition of V defined by Ea. (jl5[) . The symplectic eigenvalues v oi V (the ordinary eigenvalues of the matrix 
|2ir2F|) can be expressed in terms of the symplectic invariants deta', det/3' and det7' of V, and are given by 



V = ±^2{AV ± V (Al/)2 - 4det(F)) = ±2^(^11 T t'i3)(w22 ± W24) (65) 

where AV = deta' + det/3' — 2det7', and the matrix elements wn, W13, W22, W24 are given in Ea.([M|). The smallest 
symplectic eigenvalue of the partial transposed variance matrix V is found out to be 

= 2^(^11 - '!^13)(l^22 +V24) 

^ 1 + V jl + v 4^^(1 + 7/- 20 exp(jgy ^^ 

" h + ^ + 2^eMT^)l 2r; + (1 + r;)3 

/1 + 7? , V77(2C(l-7?)^ + 4ry(l + 7;))exp(-C/(l+^ 

+ OTvf ^''^ 

2 

where ?7 = ^ and C = ^■ 

2 

In Fig.l we plot versus the two parameters V ~ ^ and C — From the graph it is clear that the two-particle 
state is separable i>_ = 1 when po = and a ^ p. As expected, the state is entangled for all other values of these 
parameters. 



B. Single particle on a non-commutative plane 



Let us now inspect the wave function given by Eq. (j37p for any possible role of the underlying non-commutativity on 
its separability property. A single particle in two-dimensions can also be viewed (if the issue of statistics is disregarded). 
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FIG. 1: (Coloronline) The smallest symplectic eigen value V- is plotted versus the dimensionless parameters 11 — % ^^'^ C ~ Is'- 



as two particles (or modes) in one-dimension, so that the wave function (|37p can be regarded as the direct product 
of two particles on a line. On the basis of this argument, we can study the separability or inseparability of the 
two-particle state described by the effective commutative variance matrix Vq in the para above Eq. ()48p in Section 
IIB. Here our task is to find the smallest symplectic eigenvalue z>_ of the partially transposed state described by the 
transposed variance matrix (V^) obtained again from the variance matrix using the transformation (fT51) . The 
smallest symplectic eigenvalue is found out to be 1, i.e., v- = 1, which docs not depend on the non-commutative 
parameter 9. Therefore, using the PPT criterion we can conclude that the two-mode state described by the variance 
matrix (|43p is separable for any non-commutative parameter 9, since as expected, non-commutativity docs not play 
any role at the single-particle level. 



C. Two particles moving on a two-dimensional commutative plane 



Let us consider two particles moving on a two dimensional plane described by the wave functions 




(67) 



(68) 



The composite wave function for the two particles is constructed to be 



*(Pl,P2) = N2[lpl{pi) ® V'2(P2) + V'2(Pl) ® "01(^2)] 



(69) 



where ~ tttt, / „ i^i,,s is the normalization constant. 

^ 2(l+cxp(— 2a aMl) 



The variance matrix can be written in the covariant form as 




(70) 
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where Aij ,Bij ,Cij ,Dij ^Eij and Gy are 2x2 block matrices given by 
Ay = (ij Xj ® 7)5, - ({xi ig) I))^{{xj ® 

= {I ® Xi Xj)q, - {{I ® Xt))4,{{I <S) Xj))q, 
r r _ 

= iV2^[2aia,- + ^ + ^exp(-2a|ap)], (71) 
a a 

B,, = ((I^!^I±^®/))^-((£,®/))^(p^.®/)^ 

= 0, (72) 

Cy = (g) ij)* - ((xj (g) /))*((/ ig %))* = -2iV|ajaj, (73) 

Aj = (ij (gpj))* - (x, »/))*(/ ®pj)vi, = 0, (74) 
^'ij = ® /)* - (g (g /)* = (/(g PiPj)^ - (/(K)pi))*(/(g)pj)^ 

= [q!% + (a(5y - 2ajajQ!^)exp(-2Q!|a|^)], (75) 

= {pi<E)Pj))^, - (pj «i <E)Pj)^, = 2Af|a,aja^exp(-2Q;|ap) (76) 

To simplify the calculations and without any loss of generality we can take 02 = 0. The 2x2 and 4x4 submatrices 
whose determinants are invariant under symplectic transformation are as follows: 

„ _ _ Mil Bn\ ^ _ « _ M22 B22\ ^ _ f Cn ^11 

"a — Pa Id 77 I ' Pb Id I i /a — In 

V^ll -^11/2x2 \-^22 ^22^2x2 V-'^ll '^11/2x2 

/ C22 £'22 \ c f eta Ja\ X ( Oib /^^x 

"^"U- G22j2x2''""l7j /3j,.,''^"l7^ ^''^ 



The matrix elements ga, Qb, Qc, gd,'ma,n^b,''nc,''nd of the standard form of the matrix (j52p are determined by the 
symplectic invariants det(Q;a) = det(/3a) = = det(ah) = det(/3f,) = g1 = g% det(7a) = marric, det(7h) = mbmd 
and det(y) = det((5a). det((56). 

Our task is to calculate the symplectic eigenvalues of the matrix V. The eigenvalues of the matrix |2zriT^| are 
invariant under the action of symplectic transformation on the matrix V. Therefore, the symplectic eigenvalues jy^'s 
can be expressed in terms of symplectic invariant quantities detaa, deto;;,, det/3a, detf3b, det7a, det7b, detSa, detSb and 
are given by 



± V 2 A, (V) ± 2 V(A,(l/))2-4det5a, (78) 



vy = ±y2A.y(y)±2^(Ay(\/))2-4det5f, (79) 

where A^; (V) = det aa + det /3q + 2 det 7a and A^ {V) = det ab + det /S^ + 2 det 7f,. The smallest symplectic eigenvalue 
for the commutative system is given by 

v'l^min{v'^_,v^_) (80) 

where 

,c 



y^2A,(F) - 2yJ{^,{V)Y - 4det<5a, (81) 



V, 



2/^y{V)-2J{^y{V)Y ^AdetSb (82) 



After taking partial transpose of the variance matrix V^'^\ The smallest symplectic eigenvalue is given by 

^min{v^_,v'^_) (83) 

where 



- 2v/(An+Cn)(i;ii-Gn) 



1 



(1 + exp(-2aa2)) 



^/(l + exp(-2aa^))(l + (1 - Aaal)exp{-2aa\)), (84) 



,~,c 



Vy_ = 2VA22-E22 = 1 (85) 
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FIG. 2: (Coloronline) The symplectic eigenvalue is plotted versus aai. 



It is clear from Fig. 2 that the smallest symplectic eigenvalue is less than 1 except at ai = and ai oo. For 
these two cases the smallest symplectic eigenvalue takes the value 1 and hence the state ^'(p) given in Eq. ((69)) is 
separable. For all other cases the state is entangled, as is evident also from the structure of the wave function in 
Eq.(|69l). 



D. Two particles on a non-commutative plane 



Let us finally consider the two particle state 5* relevant to study the non-commutative dynamics on a two- 
dimensional plane: 



(86) 



where Fg ^ = exp (j^—^){pi 'S' Pj)^ is the (inverse) twist operator, as defined earlier in Section II, and is a normal- 
ization constant, tpi and "02 are given by 



-= exp — p -—y-ip-a ^— 

na 2a 2 4 



(87) 



and V'2 = V'l (p; - a, -Po , a) ■ 

The variance matrix can be obtained after a long but straightforward computation (see Appendix II and III) in the 
covariant form and is given by 



V 



( -^ij ^ij ^"ij "^"ij 

-Bji -^ij -^j'i ^'ij 

Cij -^ij '^'ij 

\D - - - 



ji Gij Bji Eij / 



in the basis ^ = {x'"i \ x'"2\p'i\p'2\^^i\ ^^2\Pi\p2^)^ where Ay, , Cij, Dij, Ey and Gij are 2x2 block matrices 
whose explicit forms are given in Appendix II. 

In order to investigate the presence of entanglement in the two-particle state ^, described by the variance matrix 
(|88p . we have to again use the transformation (j34p which transforms the variance matrix for the state in the non- 
commutative plane to the effective variance matrix (corresponding to the variance of commutative observables) in 
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the commutative plane. Thereafter, we can apply the partial transposition to detect the entanglement in the state ^' 
using the PPT criterion. Now, the transformed variance matrix in the commutative plane is given by 







Cij 
















Cij 


















8x8 



"^^X n\ n'. T. u'. (89) 



where Aij, Bij, Cij, Dij, Eij and Cij are 2x2 block matrices related to the earlier (barred) matrices as follows: 



Aij — Aij + - {ejkBik + tikBjk) + —eik^jiEki 



B,j - Bij + -tikEkj 



Cij = Cij + -{ejkDik + iikDkj) + -^^ik^jiGki 



(90) 
(91) 
(92) 



Dij - Dij + -e^kCkj 



Cij — Cij 



E - — E- 

'^3 — *J 



(93) 

(94) 
(95) 



We first determine the physicality condition for the variance matrix (j89p . and then perform the partial trans- 
position of it. Finally, we compute its symplectic eigenvalues. For a straightforward illustration of the effect of 
non-commutativity on entanglement, let us consider the following case. We assume that the two particles described 
by the wave functions -01 and ^2 are in a static position with respect to each other, i.e., ~ Q and that one of the 
components of h (defined as bi = ai + OijPoj = Ui) is zero, (say, 62 = 0). With the above assumptions, the symplectic 
eigenvalues 1^^^ and i^yj^' are obtained to be 



,NC 



,NC 



612 



2V(Aii-Cii)(Sii-Gii) = 2\/((Aii - Cii) + -(£^22 - G22))(Sii - Gn) 



2N^ 



V 



302^2 4(3a202+4)exp(^±^) 

— h 46ja H -^7T^ — 

4 ^ {a^9^+4)^ 



16exp( (^4"+4) ), 



-y^ . , 

^ (a26)2 + 4)3 

2V(^22-G22)(S; 



- 16a52 + 4), 



G2 



2\/(A22-G2 



-{Ell — Gii))(i?22 — G22) 



2A^2 



V 



22 4 -X f ) 
39^ ^ expl(^.g.+4)i ^3g,^4 ^ ^ _ 3202^352) X 



(0,25/2 _^ 4)3 



16exp((^^5g^) 



-4)2 



(96) 



(97) 



The variance matrix ([89|) represents the physical state if it satisfies the appropriate physicality condition for a 
non-commutative system explained in section IIC, viz. 



i.^^r = m^n((.f_^)^ «_^)^) > i.!^^\e))l 



(98) 
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In order to extract the ^-dependence in the physicaUty condition, we have to first obtain the minimum uncertainty 
state for which the function AiiAx2(^) is minimized. The above function can be read off from the elements of 
the effective commutative variance matrix (j89|) . and is given by A11A22 in terms of parameters defined in Eq. (|90p . 
Note that {0))min depends also on the state parameters a and bi. For the present computation we minimize 
AxiAx2{0) with respect to the spread a for various values of the particle separation oi ~ bi. Or in other words, 
choosing a particular state corresponding to a specific value of 61, we compute the value of a (say, amm) for which the 
uncertainty Aa;iAa;2(^) is minimized. Next, we find the minimum symplectic eigenvalue {v^^ {0))min by substituting 
the values of amin and bi in Eqs. (|55|) and ([57]) . We then repeat this computation for other values of &i. Wc find that 
the value of minimizing the position-position uncertainty stays constant for a specific range of the parameter 

61, and then again acquires a different contant value for another range of 61. In this way we obtain the physicality 
condition for ranges of states parametrized by various values of particle separation bi. The values of {1^-^ {0))min 
corresponding to three different ranges of values for the parameter 61 are plotted in Figs. 3, 4 and 5, with the respective 
minimum uncertainty conditions, i.e., a,„i„ ~ 1-4/^, amin — 1.6/0 and Umm — 1.8/0 displayed respectively on the 
figures. 




FIG. 3: (Coloronline) The smallest symplectic eigenvalue {Vx^^Y of the partially transformed variance matrix V and the function 
iy^'^ {9))min representing the entanglement criterion are plotted versus obi. Here amin — lA/9. The state is entangled since 
{Vx^Y < {&))'min for all values in the parameter space. 



Next, in order to apply the PPT entanglement criterion, we perform partial transposition of the variance matrix, 
which corresponds to time reversal, or in terms of continuous variables, sign change of momentum variables. Thus, 
when we take the partial transpose, we have to change the sign of all the components of the momentum variables of 
one of the particles. Therefore, the variance matrix V transforms under partial transposition as 

V = KVk (99) 

where A = diag{l, 1, 1, 1, 1, 1, —1, —1). After performing partial transposition operation, we compute the eigenvalues 
z>_ of the matrix |2if2y|. Due to the partial transposition operation, the symplectic invariant quantities detja and 



16 




FIG. 4: (Coloronline) The smallest symplectic eigenvalue {C'xf^)^ of the partially transformed variance matrix V and the function 
{u^^' {6))min representing the entanglement criterion are plotted versus abi. Here amin — 1.6/6. The state is entangled since 
{u^f')^ < {f-'^ {6))1nin for all values in the parameter space. 



det^b just flip their sign. Hence, the expression for the smallest symplectic eigenvalue can be written as 

z>5^ = rr^^n(z>f5,^;^^_^) (100) 

where 



- \l2A,{V) - 2^{A,{V))^ - 4detJ„ A, )= det aa + det/?a - 2 det7a, (101) 
^2Ay{V) - 2^{Ay{V))^ - 4det6t, Aj,(l/) = det + det /3b - 2 det7f, (102) 



~NC 

y- 



For a non-commutative system the state is entangled when the smallest symplectic eigenvalue = given in 
Eg. () loop satisfies the condition (following from Eq. ([59|l ) 

(^~_^C)2 < (,5C(^))2^^^ (103) 

We now perform the partial transposition of the variance matrix ([89|) and compute the symplectic eigenvalues i>^^ 
and of the partially transposed matrix with the help of Eqs. ()10ip and ()102|1 . These are given by 

on2 2 4(80*a'^b'i + Se^a'^ + WO^a^ + 16)exp( ) 



3 



4 (a26l2 + 4) 

(a202 + 4)4 Jil+ (a202+4)3 J' (104) 
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FIG. 5: (Coloronline) The smallest symplectic eigenvalue {C'xf^)^ of the partially transformed variance matrix V and the function 
{u^^' {6))min representing the entanglement criterion are plotted versus abi. Here amin — 1-8/6. The state is entangled since 
{u^f')^ < {f-'^ {6))1nin for all values in the parameter space. 



where iV^ = 2(i+4/5oxp(-8ab^/5)) • "^^^^ state is entangled if the smallest symplectic eigenvalue satisfies the condition 

(fTua . 

The smallest symplectic eigenvalue {v^^Y is plotted versus ah^ in Figs. 3, 4 and 5 respectively for different ranges 
of the parameter bi corresponding to the respective displayed values of amin- From the figures it can be seen that 
v^'^ satisfies the condition p03p . for a range of parameter values thus showing that the state is entangled for these 
parameter values. 

Now we are in a position to compare the amount of entanglement in the commutative (i.e., taking the = limit in 
the composite state of Eq. ([55)) ) and non-commutative (i.e., 6* 7^ 0) case, respectively. The respective magnitudes of 
entanglement can be quantified using the logarithmic negativity [2^. The two expressions can be formally written as 

= max[Q,-]^log2[v'l)\ (106) 
E^'^ = max[Q,-]^log2[v^f] (107) 

where = {i'^^)^ — i'^-'^ {(^))min + ^'^d the smallest symplectic eigenvalues and are obtained from 

Eqs.dHa]) and (fTOOl) . respectively. 

In Figs. 6, 7 and 8 the magnitude of entanglement E is plotted versus abi respectively for different ranges of 
bi (corresponding to the appropriate values of amin for the respective ranges of &i) for the commutative as well as 
the non-commutative case separately. One can see that the underlying non-commutativity leads to a reduction of 
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FIG. 6: (Coloronline) The logarithmic negativity E'^ and E'^'~" are plotted versus abf respectively for the commutative and 
non-commutative systems. Here Omin — lA/6. 



entanglement in all the cases. This feature is most prominent, i.e., the difference in the magnitute of entanglement 
between the commutative and the non-commutative system is maximum, when ah\ ~ 0(1). For larger separation 
of the two particles (61 = ai for the case we have studied), this difference tends to decrease, and the state becomes 
separable for asymptotically large values of bi in the commutative as well as the non-commutative case. This clear 
manifestation of non-commutativity hindering the entanglement of two particles at small distances, could presumably 
be correlated to the enhanced variances of the physical observables coming from the modified space-space and phase- 
space uncertainty relations (pS)) and (|T71) in non-commutative space. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have investigated the role of position-space noncommutativity on the entanglement of bipartite 
systems. We focus on Gaussian states in the continuous variable position-momenturn phase space, where the formalism 
for detecting and quantifying bipartite entanglement is well-established [13, E, H^. We have first extended the 
symplectic formalism for studying entanglement for a bipartite Gaussian state residing on a two- (spatial) dimensional 
plane. The validity of this formulation is illustrated with the help of an example of a two-particle state which is shown 
to exhibit expected entanglement properties. Wc then prescribe our formulation for extending the PPT criterion 
for entanglement of non-commutative systems. The ^-dependence of the entanglement criterion is derived. As an 
off-shoot wc show that for a single particle on the two-dimensional noncommutative plane, whose state may be viewed 
as a composite state of two one-dimensional modes, the underlying non-commutativity does not impact the a priori 
separability of the modes. We have finally presented an example of a non-commutative state of two particles. This 
represents a pure state and is manifestly entangled in the commutative {9 = 0) limit. As a striking feature of our 
formalism, we are able to show that the magnitude of entanglement reduces under the effect of non-commutative 
dynamics for a whole range of physically allowed parameter values. 
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FIG. 7: (Coloronline) The logarithmic negativity E'^ and E'^'~" are plotted versus abf respectively for the commutative and 
non-commutative systems. Here Omin — 1.6/9. 



A notable consequence of this study has been our result of how the criterion of entanglement itself is modified 
in the presence of non-commutativity. Our goal in this paper has been not to regard entanglement as a resource 
for performing information processing tasks, but rather to study the impact of noncommutativity on entanglement 
itself. Indeed, the role of (anti-)symmetrization associated with the statistics of indistinguishable particles on their 
entanglement is itself not yet well-understood, and any entanglement resulting from such (anti-) symmetrization might 
not be amenable for use as resource [s^. It may be relevant to observe here that quantum entanglement is difficult 
to preserve under the interactions of the system with the environment. Decoherencefs^ is an effective phenomenon 
where one is practically unable to monitor the huge number of degrees of freedom of the environment. It remains an 
open question whether the loss (or reduction) of quantum coherence could be an intrinsic feature of the dynamics of 
individual systerns, instead of being a purely effective one. There are indeed some models such as gravity induced state 
vector reduction [Sg and continuous spontaneous localization [3^, that produce quantum decoherence as an ingredient 
of the basic quantum dynamics itself, though these schemes are yet to be verified experimentally. On the other 
hand, the idea of space-space noncommutativity is well-inspired at the fundamental level, and is certainly an off-shoot 
of the physics in strong gravitational fields in the presence of event horizons of black holesQ. Non-commutative 
dynamics being able to restrict the entanglement of two particles at short distances is a remarkable consequence of 
the analysis presented in this paper. Such an effect of non-commutativity is expected to have interesting consequences 
on the entanglement of quantum states and related information capacities wherever position-space noncommutativity 
is effective. 

Acknowledgments: We would like to thank Subhash Chaturvedi and Frcderik G. Scholtz for discussions, and Soumen 
Mandal for help using Mathematica to compute the integrals for the non-commutative variance matrix elements. ASM 
would like to acknowledge support from a project funded by DST, India. 
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FIG. 8: (Coloronline) The logarithmic negativity E'^ and E'^'~" are plotted versus abf respectively for the commutative and 
non-commutative systems. Here Omin — 1.8/9. 



Appendix I: Standard form of the variance matrix 



By the following steps one can bring a variance matrix into the standard form from which its sym- 
plectic eigen values could be easily obtained. Step-I: The variance matrix can be expressed in the basis 
{x^i\ Pi^^ ,x^i\pi'\, , P2^^ ) in terms of 2 x 2 block matrices as: 



/ 



V = 



.(1) 



pi'' 



(1,2) \ 
V ' 



7 

(1,2 

Tv ■ 



(108) 



8x8 



The superscripts 1 and 2 refer to the first and second particle, respectively, and subscripts xi and X2 refer to the 
two perpendicular axes a: and y, respectively. Using Williamson's theorem [24| . we can choose two local symplectic 
transformations 81,82 G 5*^(4, i?) which transform 



V^{8i(^ 82)V{8'[ (g) 8j) = V 



( gal 





9bl 



((^,a,2)v^T 







U'^^'y 


gdi J 



(109) 



where / denotes the 2x2 identity matrix and denotes the 2x2 null matrix. As a result of the transformation, the 
cross-correlation block matrices are modified, and they are denoted with prime. 
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Step-II: There exist symplectic transformations S3, € Sp{2, R) which transform the variance matrix V as: 

/S^ 



v 



/'S3 

5*4 

53 

V 54, 



v 



5j 











(110) 



V 5j 



The transformation brings the 2x2 block matrices into diagonal form as follows: 



rria 
rrir 



Further, the transformation does not effect the 2x2 block matrices of the type 



9^ 

9^ 



(111) 

{i = a,b,c, d) since they 



are proportional to the identity matrix. 

After performing these steps, the standard form of the variance matrix can be expressed in the form 



/ 9a 











rria 





qa 


\ 





9a 
















qb 








9b 





qa 





7716 














9b 





qb 





md 


ma 










9c 














rric 





qb 





9c 








qa 
















9d 





\ 


qb 
















.9d / 



(112) 



8x8 



Appendix II: Non-Commutative variance matrix elements 

The entries of the variance matrix V for the non-commutative case are obtained from computing the following 
relevant expectation values: 



2r I XiXj -\-XjXi 



)V'2 



Q ^ X/ ' Xf ' I ■ ^ ■ ^ 

+ {<^ikPkejiPi),p2} + 2-Re(( ' ^ ^ ^ ' ^ -f)-F0"^To)^i«.V 

+ Yi?e((/ (g) eifcPfcej/p/)Fg"^fo)^i,»^2 + 9Re{{xi ® ej7Pz)-Fe~^To)v.i<»iA2 
+9Re{{xj (g) eikPk)Fg'^To)^^(g,^^], 



(113) 



X ' T) ' I • T) OC ' 

2 (fe><Ai (e»fePfe>V'2 - {^ikPk)4>i {Pj)4>2) + 2i?e(( ' ^ ^ ^ ' ® I)Fg'^To)^^^ 
+9Re({pj (g) ejiPi)-^e'^^o>Vn«>V2] 



(114) 



+2Re{{piPj (g) /)-Fe"^ro)^i«,^2] 



(115) 
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0^ X ' X ' ~\~ X ' X ' 

+^i{<^ikPkejiPi),pi + {eikPkejipi),p2) + 2i?e((/(X) ' ^ ^ ^ ' )-Fg~^To)^i^^, 

+—Re({e^kPk^jiPi ® I)Fg''^To)^^(g,^^ ~ 9Re{{eikPk a;j)^'e"^To)^i<8.V'2 

-9Re{{ejipi (g) Xi)Fg"^fo)^i@)^2] (116) 

X ' J) ' ~\~ TJ ' X ' 

2((Pj> !Ai(e»fePfe>V'2 - (e»fcj5fc)^i(pj)^J + 2Re{{I(» ' \ ' ' )Pi^^o) ,pi0i,2 

-ORe{{eikPk '^P])Fg'^TQ)^^^^^] (117) 

{Fg'\l®p^pj)Fe)^ ^ N\p.,pj)^, + (^Pi)^2 

+2Re{{I (g) Pi]5j)i^e"^Tn)^ie5V2] (118) 

{Fg'^ixi (g)Xj)^e)* = iV^[(xi)^i(lj)^2 + {^])4'i{^i)i>2 + -j^{-{^i<^jkPk)ii^x - 

{xjeikPk)i>i) + {xiejkPk)i,2 + {xjeikPk)i,2) + -^{{P^)i>l{P])i>2 - {Pj)i>APi)i>2) 
+2Re{{xi (g) :r-,)Fg~^fo)^i,»V'2 + (fRe{{I ® iikPkXj)Fg'^fQ)^^(^^^ 

-9Re{{xiejkPk <^ ^)-Fir^To) Vi«>^2 + y ^e((Pj «) Pj)-F0"^to)i>i(»^2] (119) 

9 

{Fg'^{x^(g)pj)Fe)-i> = N'^[{xi)^,AP])i'2 + {Pj)^A^i)^2 - 2{{<^ikPkPj)i,i - 

{(-ikPkP 3)^,2) + 2i?e((lj (g) Pj)Fg'^fo)^^^^^ + 9Re{{I ® eikPkP3)Fg'^fo)^^^^^] (120) 
(Fg"^(K®lj)^e),i, = N'^[{x3)^,APi)i'2 + (Pj)^i(^j)^2 - 2{{<^jkPkPt)i,t - 

{i3kPkPi)^2) + '^ReiiPi ® ^j)-^g~^fo),/;i(»V'2 - 9Re({ejkPkPt ^)-F'e"^To)^i«u/;2] (121) 
(F(r^(pj (gpj)Fe)* = A^^[(p»)vi(Pj)'A2 + (Pj>^i(Pz)v2 

+2i?e((K ® p,)/'g-2^o)v.i®^2] (122) 
Using the above relations the elements of the CM given by Eq. (j88| are obtained to be 

A3 = (F-i( (^^ + ^^) ® - {Fg-H■^^ ® /)F^^>*(i^,-\l^, ® /)F«)* 

= lllEI±lM)j^,)^ _ (F-i(/®i,)/',)^(i^^'i(/®4.)/',)^ 

2 exp( ) exp( 4^ ) K^^^^- W + 4a^6,;^, + 4a5.,) - 
49^a'^{ejkbkPQi + eikbkPoj) - W"^ {p^iPoj + 2a^eikbk>^3ibi - 3a(5y ) - 



2-0 , 8a^e% . Hpli+Ph) , .zMfcMl^W 

( — ^.-fcPOfe + (^202+4)2 ^ ) eM (^202+4) ))]' 

iV'[— Po.e.fcm - 2aey +Po,6. + 8cxp( ^"^^ cxp( (^2g2^4) ) x 

[{6l3a^(2aejfe6fc6i - e^j) + 29^a^po,bj + Aea{2abje,kbk - e,j) + Spoi^j} x 

N^i^ibje^kPok + b,e,kPok) - 26,5, + 2 exp(^^^21±Zo2) ) ^ 

'^^P^ (^g2 ^4)^^ )[i^^"'^(Po»Po3 + '^a^bibj) + Ae'^a^hki^ikPO] + fj/cPoj) 
86'^a^(poiPoj(l + 2a^) + 2a^eikejibkbi) + 169a^bk{eikP0] + ^jkPai) + 16poiPoi} 
+ 4,-3,1 - «^(^,,p„, + ^^^^ exp(d£k±l4), , 
-8a(62 + 62) -0 8a202^j ,-(Pol+Po2)^ , -8a(6? + ^^1) 

{F-\-i-®pj)Fg)^ - {F-^{Wi(E)I)Fg)^{I(g)pj)^, 

A^'I^Poj e^fcPofe - Moj - 16cxp( cxp( ^^^^^^p x 

[{~e^a%,ejkbk + ^^a^Po^bJ + AOaH.e.kbk + Ap^^b,}{{a'B^ + 4)-3}] _ 

^ (^-'^P°^ + WTIF----P( 2. (^4^ +4) )) ^ 

16a'g£j-fc6fc „.„, -(Poi +Pg2) ^ , -8«(&i + ^1) 

(a202 + 4)2 expl 2^ jexp( ^^2^2 + 4) iiJ' 

{F-\p,®pj)Fg)^ - {p, ® /)vt(/ (8) Pj)^ 

. 4,-3,, - o.p,d£l^,„.p,^||±|l,, 

^^202 + 4)2 cxpl 2a (a202 + 4) 

{piPj (g)/)* - (p, (E)I)^(pj igi/)* = {I ®PtPj)^ - {I ®Pi))m{I ®Pj)m 

+ + 16cxp( )cxp( (^2g2+4) ) X 

[{e^aH2ae,kbkejibi + S,,) - 4a(2a6,,6, - S,^)}{{a^e^ + 4)-3}] - 
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where the normahzation constant is given by 



7V~2 _ 2(1 + 1 cxp(^^^il±^)exp(^^^^^i±^)) (129) 



Appendix III: Method of the calculation of various integrals 

The method followed to calculate the terms like i?e(((;r]")^ (g) /)-Fy^fo)^i(gi^2 required to obtain the expectation 
values is as follows. The generic element looks like ((A® B)F^'^To)^^(g,^2- The above expectation value can be 
calculated by inserting the complete set of momentum eigenstates. Therefore, 



({A(g,B)Fg-^To)^,»4„ = J d'kdH exp(^%fc,^,) (fc|A|?Ai)*V2(fc)(/|B|^2>*V'i(0 

d^kd^l cxp{ik,q,) f{k)g{l) (130) 



where we define Oeijlj — qi, f{k) = {k\A\^i)*^2{k), g{l) ~ (/|i?|V'2)*V'i(0- Further, Eq. (|130p can be simplified as 

d^kdH exp{tk,q,) f{k)gll) = ^ fd^kfikMk) (131) 



^2 . 

where d^l = 4^ and g{k) is the Fourier transform oig{q) = g{l) given by 



9{k) I cxp{ik,qi) g{q) (132) 
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